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Subtitle 

Stabilization  of  an  equilibrium  point  of  a  nonlinear  system  by  linear  feedback  is  studied  in 
two  critical  cases,  corresponding  to  the  presence  of  either  a  zero  eigenvalue  or  a  pair  of  pure 
imaginary  eigenvalues  for  the  linearized  system.  In  each  case,  the  critical  eigenvalues  are 
assumed  uncontrollable,  and  stabilization  is  studied  using  bifurcation-theoretic  calculations. 

ABSTRACT 

Linear  feedback  stabilization  of  nonlinear  systems  is  studied  for  systems  whose  linearization 
at  an  equilibrium  point  possesses  a  simple  critical  mode  that  is  uncontrollable.  The  results 
complement  previous  work  on  the  synthesis  of  nonlinear  stabilizing  control  laws.  The  present 
work  addresses  continuous-time  systems  for  which  the  linearization  has  either  a  simple  zero 
eigenvalue  or  a  pair  of  simple  pure  imaginary  eigenvalues.  Both  the  stability  analysis  and 
stabilizing  control  design  employ  results  on  stability  of  bifurcations  of  parametrized  systems. 


Keywords:  Stabilizers,  nonlinear  systems,  controllability,  feedback  control. 

1  Introduction 

Feedback  stabilization  of  nonlinear  control  systems  is  a  subject  which  has  received  significant 
attention  in  the  control  literature  in  recent  years.  Local,  semiglobal,  and  global  stabilization 
problems  have  been  considered.  In  this  paper  we  address  problems  of  local  stabilization, 
using  linear  feedback,  of  an  equilibrium  point  of  a  nonlinear  system 

*  =  /(*,«)•  (1) 

Here,  x  <G  1R"  is  the  state,  u  6  IR.  is  the  control,  and  f  is  sufficiently  smooth  in  x  and  u  with 

/(0,0)=0. 

Thus,  we  consider  the  question  of  existence  of  a  linear  feedback 


u  =  u(x)  —  kx 


(2) 
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( k  €  ]R1Xn)  rendering  the  origin  of  Eq.  (1)  locally  asymptotically  stable.  In  addressing 
this  question,  we  take  an  approach  which  is  by  nature  constructive.  By  standard  theory, 
stabilizability  of  the  linearization  of  (1)  is  sufficient  for  local  stabilizability  of  the  nonlinear 
system  (1).  Thus,  in  this  paper,  as  in  many  previous  studies,  we  assume  that  the  linearization 
of  (1)  at  the  origin  has  uncontrollable  critical  eigenvalues.  An  eigenvalue  is  critical  if  it 
lies  on  the  imaginary  axis.  More  specifically,  we  consider  two  critical  cases.  In  the  first, 
the  linearization  has  an  uncontrollable  zero  eigenvalue.  In  the  second,  an  uncontrollable 
complex  conjugate  pair  of  purely  imaginary  eigenvalues  occurs.  The  recent  book  (Bacciotti, 
1992)  and  the  review  paper  (Sontag,  1990)  provide  overviews  of  many  results  on  nonlinear 
stabilization.  The  papers  (Andreini  et  a I.,  1989),  (Bacciotti  and  Boieri,  1990)  address  linear 
feedback  stabilization  of  nonlinear  systems  with  uncontrollable  critical  modes. 

Since  the  right  side  of  (1)  is  smooth,  we  may  write 

x  =  /( x,  0)  +  ugi(x)  +  u2g2(x)  +  u3g3(x)  H - .  (3) 

We  can  also  expand  the  zero-input  system  corresponding  to  (1): 

x  =  f(x, 0) 

=  Ax  +  B(x,  x )  -f  C(x,  x ,  #)  +  •••.  (4) 

Denote  the  Taylor  expansions  of  the  input  kernel  functions  gi(x),  i  =  1,2,...,  appearing  in 
Eq.  (3)  as 


gd'-r)  =  hi  +  Atx  +  Bi(x,x)  +  Ci(x,  x,x)  +  ■■  ■ .  (5) 

In  (4)  and  (5),  B(x ,  x),  C(x,  x,  x),  and  Bi(x,  x),  Ci(x ,  x ,  x),i  —  1,2,...,  are  vector  quadratic 
and  cubic  forms  induced  by  symmetric  bilinear  and  trilinear  forms,  and  the  dots  denote 
higher  order  terms  (Abed  and  Fu,  1986),  (Abed  and  Fu,  1987),  (Iooss  and  Joseph,  1990).  To 
conform  with  standard  notation  in  linear  control  theory,  we  shall  henceforth  denote  b\  by  b, 
viz., 

b-~  V 

As  noted  by  Brockett  (1983),  the  only  situations  in  which  the  smooth  local  stabilizability 
question  for  (1)  cannot  be  addressed  based  only  on  considerations  for  the  linearized  system 

x  =  Ax  +  ivu  (6) 

are  those  in  which  (6)  possesses  at  least  one  uncontrollable  critical  mode.  That  is,  the  cases  in 
which  the  determination  of  existence  of  a  stabilizing  smooth  feedback  requires  one  to  consider 
the  influence  of  the  nonlinear  terms  in  (1)  are  those  in  which  the  matrix  A  has  at  least  one 
uncontrollable  eigenvalue  with  zero  real  part.  These  are  referred  to  as  the  critical  cases 
in  nonlinear  stabilization.  This  terminology  derives  from  the  classical  nonlinear  stability 
literature,  in  which  a  critical  case  corresponds  to  the  vanishing  of  the  real  part  of  at  least 
one  system  eigenvalue  (Zubov,  1964). 

This  having  been  said,  we  wish  to  note  that  there  are  cases  in  which  a  system  is  stabiliz- 
able  by  linear  methods,  but  a  stabilizing  controller  designed  using  nonlinear  considerations 
is  preferable  to  a  traditionally  designed  linear  controller.  For  example,  a  system  may  pos¬ 
sess  eigenvalues  on  the  imaginary  axis  which  are  nearly  uncontrollable.  In  such  a,  case,  a 
linear  feedback  giving  an  adequate  margin  of  stability  may  have  an  unacceptably  large  gain, 
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whereas  there  might  exist  a  nonlinear  controller  of  low  gain  which  also  stabilizes  the  system 
adequately.  There  is  evidence  (Marino  and  Kokotovic,  1986)  that  high  gain  linear  feedback 
which  achieves  a  large  margin  of  stability  for  the  linearized  system  may  compromise  stabil¬ 
ity  of  the  nonlinear  system.  A  careful  investigation  of  these  issues  would  require  a  separate 
effort. 

Several  approaches  have  been  used  in  addressing  stabilization  problems  for  (1)  in  critical 
cases.  Following  the  introductory  work  of  Aeyels  (1985),  center  manifold  reduction  has  been 
frequently  employed  in  the  literature  on  feedback  stabilization;  see,  for  instance,  (Behtash 
and  Sastry,  1988),  (Andreini  et  ah,  1989),  (Bacciotti  and  Boieri,  1990),  and  (Boothby  and 
Marino,  1989).  An  alternative  but  mathematically  equivalent  approach  was  introduced  by 
the  authors  in  (Abed  and  Fu,  1986),  (Abed  and  Fu,  1987).  Here,  a  bifurcation-theoretic 
framework  is  employed,  and  stabilization  of  an  equilibrium  point  in  a  critical  case  is  tied 
to  the  stabilization  of  bifurcated  solution  branches.  This  approach  is  also  employed  in  the 
present  paper.  This  approach  affords  the  advantage  of  calculations  performed  directly  on 
the  original  system  dynamic  equations. 

From  a  control-theoretic  viewpoint,  it  is  interesting  to  consider  the  feasibility  of  employing 
the  common  linear  feedback  laws  in  the  local  feedback  stabilization  of  nonlinear  systems. 
This  work  is  also  motivated  by  Theorem  3  of  Abed  and  Fu  (1987),  which  asserts  the  following: 
Under  certain  generic  conditions  (see  (Abed  and  Fu,  1987)),  the  local  feedback  stabilization 
problem  for  (1)  whose  Jacobian  A  has  an  uncontrollable  zero  eigenvalue  is  not  solvable  by  a 
smooth  feedback  control  with  vanishing  linear  part. 

In  essence,  the  purely  nonlinear  stabilizing  feedback  laws  in  (Abed  and  Fu,  1986),  (Abed 
and  Fu,  1987)  are  constructed  by  ensuring  negativity  of  certain  so-called  bifurcation  stability 
coefficients,  the  values  of  which  can  be  determined  using  formulae  given  directly  in  terms  of 
the  system  dynamic  equations.  It  appears  that  these  bifurcation  stability  coefficients  play  a 
role  similar  to  that  of  dominant  open-loop  pole(s)  in  linear  control  theory.  Thus  the  problem 
investigated  here  may  be  viewed  as  a  corresponding  “nonlinear  pole”  assignability  problem. 

The  strategy  for  studying  stabilizability  employed  in  previous  papers  of  the  authors  (e.g., 
(Abed  and  Fu,  1986),  (Abed  and  Fu,  1987),  (Fu,  1990),  (Fu  and  Abed,  1993))  will  also  be 
used  here  to  study  existence  and  synthesis  of  linear  stabilizing  feedback  controllers.  The 
analysis  is  more  involved  than  that  in  (Abed  and  Fu,  1986),  (Abed  and  Fu,  1987),  however. 
This  is  because  linear  feedback  modifies  the  Jacobian  matrix  A,  unlike  the  case  with  purely 
nonlinear  feedback.  This  results  in  two  important  differences  between  the  work  in  (Abed 
and  Fu,  1986),  (Abed  and  Fu,  1987)  and  the  present  endeavor:  First,  the  linear  feedback 
gain  vector  k  must  now  be  restricted  t.o  ensure  that  the  nominally  stable  eigenvalues  remain 
stable  for  the  closed-loop  system.  Second,  in  (Abed  and  Fu,  1986),  (Abed  and  Fu,  1987)  the 
pertinent  closed-loop  bifurcation  stability  coefficients  are  polynomial  functions  of  the  control 
gain.  In  contrast,  we  shall  be  faced  with  rational  functions  of  the  gain  k  in  the  present  paper. 

If  it  happens  that  b  =  0,  which  occurs  in  particular  if  (1)  is  affine  in  the  control,  then  A 
is  invariant  under  any  linear  feedback.  Under  this  condition,  the  complications  just  noted 
do  not  occur.  The  simplified  results  for  this  situation  follow  easily  from  the  more  general 
results  presented  here. 

The  remainder  of  this  paper  proceeds  as  follows.  In  Section  2  the  basic  problem  set¬ 
up  is  given,  including  expressions  for  the  stability  coefficients  of  the  open-  and  closed-loop 
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systems.  In  Section  3  the  stability  of  the  linearized  closed-loop  system  is  considered.  Section 
3  also  addresses  the  computation  of  some  relevant  matrix  inverses  as  well  as  of  the  closed- 
loop  eigenvectors  associated  with  the  critical  eigenvalues.  The  main  results  of  the  paper  are 
presented  in  Sections  4  and  5:  Section  4  contains  the  results  in  the  case  of  an  uncontrollable 
zero  eigenvalue,  while  Section  5  addresses  critical  systems  with  an  uncontrollable  pair  of  pure 
imaginary  eigenvalues.  Concluding  remarks  are  collected  in  Section  6. 

Notation.  In  what  follows,  ||-||  stands  for  the  Euclidean  norm:  ||w||  =  (vTv)i  for  v  £  JR"  and 
||i>||  =  ( vHv)*  for  v  £  Cn.  The  real  and  imaginary  parts  of  a  scalar  expression  are  denoted 
by  Re{-}  and  Im{-},  respectively.  A  superscript  T  (resp.,  H)  indicates  the  transpose  (resp., 
Hermitian  transpose)  of  a  matrix  or  vector.  With  k  £  IR1Xn,  denote  by  7 Z-k(F(-))  the  range 
of  a  function  of  k  over  a  domain  ||fc||  <  k: 

Kt(F(-))  :=  {F(k) :  ||*||  <  k).  (7) 

2  Stability  Criteria  and  Stabilization  Strategy 

As  mentioned  above,  this  paper  addresses  the  stabilization  by  linear  feedback  of  the  origin 
of  Eq.  (1)  in  two  critical  cases.  In  this  preliminary  section,  we  give  the  basic  problem 
set-up  and  recall  formulae  for  quantities  known  as  stability  coefficients.  The  use  of  the 
stability  coefficients  in  stability  assessment  for  the  origin  in  the  critical  cases  studied  here 
is  summarized.  The  effect  of  linear  feedback  on  the  values  of  the  stability  coefficients  is 
determined.  Finally,  we  outline  an  overall  strategy  for  the  investigation  of  stabilizability  by 
linear  feedback. 

Consider  the  zero-input  system  (4).  In  this  paper,  the  Jacobian  A  of  this  system  is 
assumed  to  satisfy  either  of  the  following  two  hypotheses: 

(S)  A  has  an  eigenvalue  Ai  =  0,  with  the  remaining  eigenvalues  A2, . . . ,  A„  in  the  open  left 
half  complex  plane; 

(H)  A  has  a  pair  of  simple,  pure  imaginary  eigenvalues  Ai  =  iuc  and  A2  =  —  i(Jc,  with  the 
remaining  eigenvalues  A3, . . . ,  An  in  the  open  left  half  complex  plane. 

For  ease  of  reference,  say  that  Case  (S)  (resp.,  Case  (H))  prevails  for  Eq.  (1)  if  A  satisfies 
hypothesis  (S)  (resp.,  (H)). 

It  is  convenient  to  denote  by  Ac  the  critical  eigenvalue  in  either  of  the  critical  cases, 
(S)  or  (H).  That  is,  Ac  :=  0  in  Case  (S)  and  Ac  :=  iuc  in  Case  (H).  The  row  and  column 
vectors  l  and  r  are  the  left  and  right  eigenvectors,  respectively,  of  the  Jacobian  matrix  A 
corresponding  to  the  critical  eigenvalue  Ac.  For  definiteness,  set  the  first  component  of  /  to 
1  in  Case  (S)  and  to  ^  in  Case  (H),  and  further  normalize  r  to  satisfy  1 

Ir  —  1.  (8) 

1  If,  may  be  necessary  to  first,  interchange  the  first  component  of  the  state  vector  with  another  component, 
to  ensure  the  nonvanishing  of  the  first  component  of  l.  This  normalization  turns  out  to  be  more  convenient 
here  than  that  used  in  (Abed  arid  Fu,  1986),  (Abed  and  Fu,  1987),  as  will  be  seen  in  Section  3.2.  In  the 
cited  references,  the  first  component  of  r  was  set  to  unity  and  then  l  was  chosen  to  satisfy  Ir  =  1  both  in 
Case  (S)  and  Case  (H).  This  modification  does  not  affect  the  form  of  the  stability  criteria. 
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In  Case  (S),  two  stability  coefficients  /?ls,  0 2s  will  be  employed  in  the  analysis.  These 
coefficients  and  their  connection  to  the  stability  of  the  origin  of  (4)  are  identified  next. 

Fact  1.  ((Abed  and  Fu,  1987),  (Iooss  and  Joseph,  1990))  In  Case  (S),  denote 

0is  lB(r,r),  (9) 

02s  :=  l[-2B(r,  {AT A  +  lTl)~1ATB(r ,  r))  +  C(r,  r,  r)].  (10) 

Then  the  origin  of  (4)  is  locally  asymptotically  stable  (resp.,  unstable)  if  0\s  =  0  and  02 s  <  0 
(resp.,  0is  ^  0,  or  0\s  =  0  and  02s  >  0).  ■ 

Two  stability  coefficients  were  introduced  above  for  the  assessment  of  stability  in  Case 
(S).  In  the  stability  analysis  of  the  origin  of  (4)  in  Case  (H),  a  single  stability  coefficient, 
denoted  02h,  occurs.  A  formula  for  this  coefficient,  and  its  connection  to  the  stability  of  the 
origin  of  (4)  are  recalled  next. 

Fact  2.  ((Abed  and  Fu,  1986),  (Howard,  1979))  In  Case  (H),  denote 

A„  :=  2Re{/[— 7?(r,  A~1B{r, f ))  -  —  2iu>cJ)~1B(r,r))  +  jC(r,r,r)]}.  (11) 

Then  the  origin  of  (4)  is  locally  asymptotically  stable  (resp.,  unstable)  if  02h  <  0  (resp., 
02h  >  0).  ■ 

We  note  in  passing  that  the  quantities  0is,  02s  and  02h  can  be  interpreted  as  coeffi¬ 
cients  in  asymptotic  expansions  of  dominant  Liapunov  exponents  for  bifurcated  solutions 
of  parametrized  embeddings  of  Eq.  (4).  Various  derivations  of  these  coefficients  have  been 
given  in  the  literature.  These  employ,  for  instance,  the  projection  method  (Abed  and  Fu, 
1987),  (Iooss  and  Joseph,  1990),  harmonic  balance  (Howard,  1979),  and  Liapunov  functions 
(Fu  and  Abed,  1993).  Note  also  that  in  the  nongeneric  situations  in  which  the  coefficients 
of  Eqs.  (9)-(ll)  vanish,  typically  stability  can  still  be  assessed  by  examination  of  further 
(higher  order)  coefficients  in  the  appropriate  asymptotic  expansions. 

Due  to  the  foregoing  considerations,  assignability  of  bifurcation  stability  coefficients  by 
feedback  plays  a  role  in  critical  nonlinear  systems  (1)  analogous  to  that  played  by  pole 
assignability  in  linear  time-invariant  systems. 

Under  either  Case  (S)  or  Case  (H),  the  basic  assumptions  we  make  on  (1)  are  that:  (i)  the 
origin  is  unstable  for  the  zero-input  system  (4);  and  (ii)  the  critical  mode  is  uncontrollable 
for  the  linearized  system  (6).  Indeed,  the  results  of  this  paper  can  easily  be  re-construed  to 
apply  to  the  case  in  which  the  origin  is  locally  asymptotically  stable  for  (4)  or  that  in  which 
the  stability  of  the  origin  is  unknown.  In  these  situations,  assignability  of  the  pertinent 
stability  coefficients  takes  the  place  of  the  pure  stabilizability  question.  Achieving  ci  certain 
degree  of  stability  or  simply  guaranteeing  stability  despite  inability  to  determine  the  stability 
properties  of  the  open-loop  system  become  the  issues  of  concern. 

Assumption  (ii)  above  is  equivalent,  by  the  well  known  Popov-Belevitrh-Hautus  eigen¬ 
vector  test  for  controllability  (Kailath,  1980),  to 

lb  =  0. 
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(12) 


Generically,  the  local  stability  or  instability  of  the  origin  of  Eq.  (4)  in  Cases  (S)  and  (H) 
is  determined  by  terms  in  the  vector  field  up  to  cubic  order.  This  follows  from  Facts  1  and  2 
above.  Thus,  next  we  record  how  these  terms  in  the  vector  field  are  affected  by  linear  state 
feedback  u  —  kx. 

Substituting  the  linear  feedback  (2)  into  (1),  and  using  an  asterisk  to  indicate  quantities 


after  feedback,  we  obtain  the  closed-loop  system 

x  =  A* x  +  B*(x,  x)  +  C*(x,  x,  x)  +  •  •  •  (13) 

where  the  linear,  quadratic  and  cubic  terms  in  the  closed-loop  vector  field  are  given  by 

A*  =  A  +  bk ,  (14) 

B*(x,x)  =  5(x,  x)  +  [kx)A\x  +  b2(kx)2,  (15) 

C*(x,x,x )  =  C(x,  x,  x)  -f  (lcx)Bi(x,  x)  +  (kx)2 A2X  +  (kx)3bz,  (16) 


respectively. 

Note  from  Eqs.  (14)-(16)  that  linear  state  feedback  affects  terms  of  every  order  in  (1).  In 
particular,  with  b  ^  0  such  a  feedback  in  general  affects  the  system  linearization:  A*  ^  A. 
This  implies  that  a  large  gain  linear  feedback  may  destabilize  modes  that  are  stable  for  the 
open-loop  system.  Hence,  in  the  stabilization  of  a  critical  system  using  linear  feedback,  one 
must  ensure  that  all  eigenvalues  of  A*,  besides  the  immobile  critical  mode,  are  not  moved 
out  of  the  open  left  half  complex  plane. 

Thus,  to  stabilize  (1)  using  linear  feedback,  it  suffices  to  find  a  vector  k  that,  firstly,  does 
not  destabilize  the  nominally  stable  system  eigenvalues,  and  secondly,  ensures  satisfaction 
of  the  stability  criteria  of  Facts  1  and  2  above  for  the  closed-loop  system.  The  first  of  these 
requirements  is  equivalent  to  asking  that  k  preserves  the  original  criticality  hypothesis  (S) 
or  (H).  The  second  requirement  we  impose  is  that,  for  the  closed-loop  system,  j3^  =  0  and 
(3 2S  <  0  in  Case  (S);  and  (32h  <  0  in  Case  (H).  Here,  by  the  formulae  in  Facts  1  and  2,  these 
stability  coefficients  for  the  closed-loop  system  are  given  in  Case  (S)  by 

(3\s  —  l*B*(r*,r*),  (17a) 

P*2s  =  /*[— 2  B*(r*,{A*T  A*  +  rTr)~lA*TB*(r*,r*))  +  C*(r*,r*,r*)\,  (17b) 

and  in  Case  (H)  by 

2Re{r[-B*(r-,r*_1B*(r",f*))  -  )«'(>■*,  (,1*  -  2tu ,J)~' B'(r‘ ,r‘)) 

+  jC*(r',r*,f-)]}.  (18) 

Here,  as  in  Eqs.  (13)-(16),  an  asterisk  indicates  a  closed-loop  quantity. 


3  The  Closed-Loop  Linearized  System 

In  this  section  we  present  calculations  relevant  to  the  closed-loop  linearized  system.  The 
concept  of  restricted  matrix  inverse  is  recalled  and  used  in  the  calculation  of  the  closed-loop 
critical  eigenvectors  l*  and  r* ,  in  addition  to  other  useful  closed-loop  matrix  quantities. 
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3.1  Restricted  Matrix  Inverses 


In  Case  (S),  we  deal  with  a  singular  Jacobian  matrix  A.  This  matrix  has  a  (geometrically  and 
algebraically)  simple  zero  eigenvalue.  In  Case  (H),  the  Jacobian  matrix  is  not  singular,  but 
we  will  none  the  less  find  it  convenient  to  solve  certain  associated  linear  algebraic  equations 
with  coefficient  matrix  A  —  iiocI ,  which  is  singular  with  a  simple  zero  eigenvalue.  In  this 
subsection,  therefore,  we  consider  linear  algebraic  equations 

Mx  =  y  (19) 

where  M  £  C"xn  possesses  a  simple  zero  eigenvalue  and  y  £  <Dn.  The  conclusions  for  M  and 
y  real  in  Eq.  (19)  follow  as  a  special  case.  The  notation  l,  r  for  the  left  and  right  eigenvectors 
of  M,  respectively,  is  as  employed  in  the  foregoing  section.  The  particular  normalization  used 
in  defining  these  eigenvectors  does  not  affect  the  restricted  inverses  obtained. 

Denote  by  J\f(M)  the  null  space  of  M,  a  subspace  of  <Dn.  Denote  by  A /*(M)X  the  orthogo¬ 
nal  complement  of  Af(M)  (a  subspace  of  <Dn).  Using  the  Fredholm  Alternative,  we  find  that, 
if  both  the  domain  and  range  of  M  are  restricted  to  A f(M)x,  there  results  a  unique  solution 
to  the  equation  Mx  =  y.  Restriction  to  A f(M)1-  is  tantamount  to  imposing  the  conditions 
lx  ■=  ly  =  0.  The  corresponding  solution  operator  is  found  to  be  given  by  the  following 
restricted  inverse  on 

=  (MH M  T  Mh  (20) 

To  illustrate  the  relevance  of  (20)  in  this  work,  note  that  formula  (9),  which  gives  /32s  in 
Case  (S),  is  rewritten  using  (20)  in  the  following  compact  form: 

[32s  =  l[-2B(r,A~B(rVr))  +  C(r,r,r)].  (21) 

3.2  Closed-Loop  Inverses  and  Critical  Eigenvectors 

Observe  that,  since  lb  =  0  (the  uncontrollability  assumption),  it  follows  from  (14)  that 
l A*  =  IA.  Taken  along  with  the  normalization  of  Section  2  for  the  eigenvectors,  we  thus 
have  that  l*  —  l  in  both  Cases  (S)  and  (H).  To  compute  r*,  write  r*  =  r  +  r  and  solve  for  f 
subject  to  the  normalization  l*r*  =  1  and  to  the  associated  eigenvector  equation  A*r*  =  Ayr*. 
(Recall  that  Ac  is  either  0  or  iuc.)  Since  l*  =  /,  l*r*  =  1  implies  that  r  satisfies  If  —  0.  It 
follows  that  the  corresponding  eigenvalue/eigenvector  equations  in  Case  (S)  and  in  Case  (H) 
lead  to  A*r  =  —{kr)b  and  ( A *  —  iu>cI)f  =  — ( kr)b ,  respectively. 

Next  we  consider  solving  the  ecjuations  at  the  end  of  the  last  paragraph  for  f.  Since 
lb  =  0  and  Ir  =  0,  the  Fredholm  Alternative  applies.  Using  also  the  fact  that  l*  =  /,  the 
closed-loop  system  restricted  matrix  inverses  in  Case  (S)  and  Case  (H)  are  found  to  be  given 
by 


A*~ 

(A*  -  iucI)~ 


A~  - 


1 


1  +  kA~b 
{A  -  iujcI)~ 


A  bkA  , 
1 


1  +  k{A  —  iu)cI)~b 


—  (A  —  iujcI)  bk(A  —  iuciy 


(22) 

(23) 
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respectively,  provided  that  1  +  kA~b  ^  0  and  1  +  k(A  —  iujcI)~b  ^  0.  These  formulae  are 
easily  verified.  Their  statement  was  initially  motivated  by  the  Sherman-Morrison  Formula 
((Golub  and  Van  Loan,  1983),  Eq.  (2.1.4)). 


Motivated  by  the  Eqs.  (22)  and  (23),  introduce  the  following  class  of  rational  functions: 

f{k]p,q)  :=  ,  k  £  R1 

1  +  kq 

It  follows  that  the  closed-loop  critical  eigenvector 

r  —  S(k)A~b 
r  —  7{(k)(A  —  iujcI)~b 

where 

S(k)  :=  A~b)  =  — 

and 

7i(k)  :=  tF(k]  r,  (A  —  iuicI)~b )  =  — 

It  is  easily  verified  that  hr*  =  S(k)  in  Case  (S)  and  hr*  =  H(k)  in  Case  (H).  These  obser¬ 
vations  will  be  useful  in  Appendix  A  in  the  computation  of  fi*s  and  /?£  ,  respectively. 

In  Case  (S),  by  (25)  we  have  that  r*  is  an  afhne  vector  function  of  the  real  fraction  S(k). 
However,  in  the  next  section  the  necessary  condition  for  stabilization  /3fs  =  0  (cf.  Fact  1) 
will  be  imposed,  resulting  in  S(k)  being  constrained  to  take  one  of  only  three  possible  values. 
Hence,  r*  is  similarly  constrained  by  the  condition  j3{s  =  0.  In  contrast,  in  Case  (H),  no 
such  equality  constraint  is  imposed.  Hence,  from  Eq.  (25),  in  Case  (H),  r*  is  an  affine  vector 
function  in  the  complex  fraction  H(k ),  a  free  variable. 


xn,  p,q  G  Cnxl. 

(24) 

r*  of  A*  is  given  by2 

in  Case  (S) 
in  Case  (H) 

(25) 

kr 

,  G  IR 

F  kA~b 

(26) 

kr 

(27) 

e  c 

F  k  A  -  iuJ)~b 

Several  matrix  inverses  associated  with  the  closed-loop  linearized  system  will  be  used  in 
the  calculations  in  Case  (II).  One  is  the  restricted  inverse  (A*  -  icocI)~  of  Eq.  (23).  However, 
in  applying  Eq.  (18)  expressions  for  the  actual  inverses  A*-1  and  (A*  -  2iucI)~1  will  be 
needed.  Since  both  A  and  A  —  2 iucI  are  invertible  in  Case  (H),  direct  application  of  the 
Sherman-Morrison  formula  ((Golub  and  Van  Loan,  1983),  Eq.  (2.1.4))  yields 

A-'bkA-1 


A* 


(A*  -  2iucI) 


=  A  - 


-i 


(A  —  2itocI) 


1  +  &A-HF 

(A 


2 iujcI)  lhk{A  —  2iuciy 


1  +  k(A  —  2iu>cI')~1b 
respectively,  provided  that  1  +  kA~xb  ^  0  and  1  +  k(A  -  2iucI)~1b  ±  0. 


(28) 

(29) 


Finally,  we  introduce  the  following  notation,  which  will  be  used  in  the  statement  of  results 
on  stabilization  by  linear  feedback  in  Case  (S)  and  Case  (H).  Denote  by  k  the  largest  scalar 
for  which  the  matrix  A  +  bk  retains  the  same  number  of  stable  eigenvalues  as  the  matrix 
A,  for  all  k  with  ||&||  <  k.  Thus,  for  all  k  with  norm  less  than  k ,  the  linear  state  feedback 
u  =  kx  does  not  destabilize  any  of  the  open-loop  stable  eigenvalues. 

2  Note  that  if  l  and  r  were  specified  in  the  way  that  one  first  normalizes  [r]i  =  1  and  then  normalizes  l 
to  Ir  =  1,  then  since  (25)  would  not  guarantee  that  [r*L  =  1,  it  would  have  been  necessary  to  take  the  new 
normalizing  factor  into  account  for  applying  (17)  and  (18).  The  present  way  of  specifying  l  and  r  makes  use 
of  the  result  l*  =  l  which  guarantees  [/*h  =  1  in  Case  (S)  and  [l*]x  =  ±  in  Case  (H),  and  thus  does  away 
with  the  need  for  renormalization. 
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4  Linear  Feedback  Stabilization  in  the  Case  of 
One  Zero  Eigenvalue 

This  section  contains  the  main  results  of  this  paper  for  Case  (S);  hypothesis  (S)  is  in  force 
throughout  the  section.  The  linear  stabilizability  of  the  origin  of  (1)  is  addressed  first, 
followed  by  a  procedure  by  which  a  stabilizing  linear  feedback  can  be  constructed  explicitly. 
In  studying  the  linear  stabilizability  problem,  it  is  natural  to  assume  that  the  origin  of  the 
open-loop  system  does  not  satisfy  the  conditions  for  asymptotic  stability  /3ls  =  0,  /?2s.  <  0  of 
Fact  1. 

In  the  remainder  of  this  section,  the  objective  is  to  find  a  linear  feedback  u  =  kx  ensuring 
that  /3^s  =  0  and  <  0,  without  moving  any  of  the  open-loop  stable  eigenvalues  of  the 
linearization  out  of  the  open  left  half  complex  plane.  The  last  requirement  can  also  be 
expressed  by  asking  that  matrix  A*  is  itself  consistent  with  hypothesis  (S).  Derivations  of 
the  expressions  for  the  closed-loop  bifurcation  stability  coefficients  and  in  terms  of 
the  gain  vector  k  are  detailed  in  Appendix  A. 

4.1  Necessary  Conditions  for  Stabilizability 

From  Fact  1,  any  smooth  feedback  which  stabilizes  the  origin  of  (I)  must  result  in  a  vanishing 
closed-loop  stability  coefficient  /3^s.  Such  a  feedback  must  satisfy  further  conditions  to  be 
guaranteed  locally  stabilizing.  Here,  we  shall  ask  that  it  also  achieve  /3%  <  0,  and  that  it 
not  result  in  destabilization  of  the  open-loop  stable  eigenvalues. 

In  this  subsection,  we  determine  necessary  conditions  for  linear  stabilizability  b}f  consid¬ 
ering  the  possibility  of  achieving  fils  =  0  through  a  linear  feedback  u  =  kx,  and  then  proceed 
to  consider  the  possibility  of  achieving,  in  addition,  the  condition  <  0. 

Substituting  x  :=  r*  from  (25)  (for  Case  (S))  in  the  formula  (15)  for  B*(x,x),  and 
grouping  terms  of  like  order  in  S(k)  (introduced  in  (26)),  we  obtain 

B*(r*,r*)  =  B(r,  r)  +  S(k)Vl  +  S2(k)v 2,  (30) 

where  the  tq  and  v2  are  real  coefficient  vectors  determined  in  Appendix  A  (see  Eqs.  (A. 3) 
and  (A. 4)).  Hence,  using  (17a),  the  necessary  condition  for  stabilization  =  0  now  takes 
the  form 

0  =  ft8(S) 

=  Pis  +  S{k)(lVl)  +  S\k){lv2).  (31) 

From  Fact  1  and  Eq.  (17a),  we  see  that  (1)  is  not  linearly  stabilizable  if  there  is  no  real 
solution  S{k)  to  Eq.  (31).  Conversely,  for  any  gain  vector  k  which  linearly  stabilizes  (1),  the 
resulting  S  =  S(k)  given  by  (26)  must  solve  (31).  When  (31)  is  truly  a  quadratic  equation, 
i.e.,  when  lv2  ^  0,  we  denote  by  A  the  discriminant  of  (31): 

A  :=(/iq)2-4/?ls(/i;2).  (32) 


Clearly,  if  A  <  0,  then  (31)  has  no  real  solutions,  and  Eq.  (1)  is  not  linearly  stabilizable. 
If  A  is  defined  (i.e.,  lv2  ^  0)  and  A  >  0,  then  (31)  has  precisely  two  real  solutions.  If 


(i 


A  is  not  defined  (i.e.,  lv2  =  0)  and  if  lv\  ^  0,  then  there  is  only  one  solution,  namely 
S(k )  =  —(3is/(lvi).  If  both  lv2  and  Iv i  vanish,  then  (31)  has  either  no  solution  or  else 
any  value  for  S(k)  is  a  solution,  depending  on  whether  or  not  fts  /  0.  Disregarding  the 
degenerate  case  in  which  all  of  lv1,  lv2  and  fts  vanish,  we  have  that  the  possible  real  solutions 
of  (31)  are  any  of  <Si,  S2,  S3  defined  as  follows: 


.  „  A1/2 

15  2  2 (lv2)  ’ 


(33) 


(34) 

Below,  we  will  need  to  consider  the  special  cases  in  which  <Si,c>2  apply  and  either  A  =  0  or 
fts  =  0,  and  the  case  in  which  S3  applies  and  fts  =  0. 


The  preceding  remarks  clearly  delineate  certain  necessary  conditions  for  linear  stabiliz- 
ability  of  (1)  under  hypothesis  (S).  The  next  theorem  collects  several  of  these.  Note  that 
these  conditions  are  also  sufficient  conditions  for  existence  of  a  linear  feedback  rendering 


Pis  =  °- 


Theorem  1.  (Necessary  Condition  for  Linear  Stabilizability)  Let  (1)  satisfy  hypothesis  (S). 
Suppose  the  origin  of  (1)  is  stabilizable  by  a  linear  feedback  u  =  kx.  Then  either  fts  =  0  or 
there  is  a  linear  feedback  u  —  kx  rendering  ft  =  0.  More  precisely,  either  fts  =  0,  or  one 
of  the  following  cases  holds: 


(a)  Iv 2  ^  0  and  A  >  0,  or 


(b)  lv2  =  0  and  lv\  ^  0. 


Theorem  1  is  a  constructive  result,  when  taken  along  with  the  derivation  preceding  it. 
For  example,  we  see  that  in  case  fts  —  0,  any  stabilizing  linear  feedback  gain  vector  k  must 
either  satisfy  S(k)  =  0  or  S(k )  =  — (/u1)/(/n2)  (the  latter  being  a  possibility  only  if  lv2  ^  0). 


4.2  Sufficient  Conditions  for  Stabilizability 

Suppose  that  (1)  satisfies  the  necessary  conditions  given  in  Theorem  1,  and  let  the  feedback 
gain  vector  k  be  such  that  S  —  S(k)  is  one  of  the  Si,  i  =  1,2,3  of  Eqs.  (33),  (34).  Thus, 
having  ensured  ftq  =  0,  the  formula  (17b)  for  ft  applies.  The  detailed  calculations,  given 
in  Appendix  A,  result  in  the  following  expression: 


In  this  expression, 


PI  =  fts(k;S) 


=  Ps(S)  +  d{S)<f>(k-,S). 

(35) 

Ps(S) 

:=  {Iwq)  +  (lw\)S  +  ( lw2)S 2  +  ( hi>3)S 3, 

(36) 

d(S) 

:=  (lvi)  +  2(/u2)<5, 

(37) 

<f>(k-,S) 

:=  f(k;  A~B*(r*,r*),A~b), 

=  T{k]  A~{B(r,  r)  +  .5m  +  S2v2),  A~b ) 

(38) 
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where  the  coefficient  vectors  ivq ,  w2,  w2 ,  w3  appearing  in  (36)  are  specified  in  Eqs.  (A.5)-(A.8) 
of  Appendix  A.  (Recall  that  v\  and  v2  are  defined  in  Eqs.  (A. 3)  and  (A. 4)  of  Appendix  A.) 
Thus  Ps(S)  is  a  cubic  polynomial  in  <S,  d(S)  is  affine  in  S ,  and  <j>(k-,S)  is  rational  in  the 
elements  of  k.  The  real  fractional  map  T  was  defined  in  (24). 

Note  that  we  have  used  the  notation  (32s(k\S )  and  (f>(k]S),  even  though  S  depends  on 
k  (see  Eq.  (26)).  This  notation  is  useful,  however,  since  it  is  a  reminder  of  the  two-stage 
nature  of  the  control  design  procedure  of  this  paper.  First,  one  determines,  if  possible,  the 
values  of  S  which  would  ensure  ft^s  =  0.  Since  any  particular  value  of  S  is  to  be  achieved 
by  appropriate  choice  of  k,  this  places  a  restriction  on  the  allowable  feedback  gain  vectors  k. 
The  next  step  is  to  find,  if  possible,  a  gain  vector  k  from  among  those  achieving  S{k )  =  S, 
which  also  results  in  a  negative  value  for  ^  ,  and  which  does  not  at  the  same  time  destabilize 
any  of  the  open-loop  stable  eigenvalues. 

In  what  follows,  we  focus  on  the  achievable  values  of  fl2c.(k]  S)  where  k  is  restricted  such 
that  S(k)  is  one  of  the  Si,  i  ==  1,2,3,  given  above.  That  is,  we  study  the  assignability  of  /?£ 
by  linear  feedback  laws  achieving  =  0. 

Consider  the  goal  of  finding  a  gain  vector  k  for  which  fi{s  =  0  and  fi2s  <  0.  In  the 
foregoing,  we  have  determined  that  this  goal  is  closely  related  to  that  of  finding  a  gain  vector 
k  which  results  in  certain  desired  values  for  S(k)  and  <j>(k;S ).  The  latter  are  fractional  maps 
of  the  form  (24).  We  are  thus  led  to  consider  the  simultaneous  solvability  of  equations 
involving  such  fractional  maps.  This  is  addressed  in  detail  in  Appendix  B,  and  a  specific 
result  is  given  in  Lemma  B.l.  From  Lemma  B.l  and  the  expressions  above  for  /?*  ,  /3 £  ,  the 
next  result  follows. 

Theorem  2.  (Stability  Coefficient  Assignability)  Let  (1)  satisfy  hypothesis  (S)  and  suppose 
lb  =  0.  Fix  a  S  £  1R  for  which  flis{S)  =  0  (see  Eq.  (31)).  The  function  (32s{k\  S)  defined  on 
the  domain  of  those  vectors  k  £  1R1X"  for  which  S(k)  =  S  is  onto  if  either  of  the  following 
holds: 

(a)  Iv 2  ^  0  and  A  >  0, 

(b)  lv2  =  0  and  lv\  0,  or 

(c)  flis  =  Iv i  =  lv2  =  0  and  either  lw3  ^  0,  or  Iw 3  =  hv2  =  0  and  lu>\  ^  0.  ■ 


Note  that  it  is  also  straightforward  to  address  in  the  same  fashion  other  special  cases. 
However,  to  limit  the  complexity  of  the  presentation,  we  prefer  not  to  discuss  these  explicitly. 


As  observed  previously,  sufficient  conditions  such  as  those  given  in  Theorem  2  do  not 
alone  imply  local  st.abilizability  of  the  origin  of  Eq.  (1)  by  linear  state  feedback.  Besides 
stability  coefficient  assignability,  a  linear  feedback  also  must  not  destabilize  open-loop  stable 
eigenvalues.  The  next  result,  addresses  linear  stabilizability,  and  involves  consideration  of 
the  issue  of  nondestabilization  of  the  open-loop  stable  modes.  To  characterize  the  linear 
stabilizability  in  a  more  explicit  fashion,  we  first  note  Eq.  (35)  and  the  notation  (38). 
Theorem  3  below  results  immediately.  The  theorem  statement  uses  the  notation  k  introduced 
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at  the  end  of  Section  3,  and  the  notation  7 Z~k(F(-))  introduced  at  the  end  of  Section  1.  This 
is  in  order  to  take  into  account  the  stable  mode  nondestabilization  requirement. 

Theorem  3.  (Sufficient  Condition  for  Linear  Stabilizability)  Let  (1)  satisfy  hypothesis  (S) 
and  suppose  lb  =  0.  Then  the  origin  of  Eq.  (1)  is  stabilizable  by  linear  state  feedback  if 
there  is  a  k  with  ||&||  <  k  that  satisfies  kA~b  ^  —1  and  S(k)  =  Si  €  '1Z^(S(-))  (defined  in 
Eqs.  (33),  (34))  for  some  i  €  {1,2,3};  and  for  this  k  and  this  St,  <f>(k]  Sj)  €  TZ‘k((f>(-,  Si)) 
satisfies  the  inequality 


Ps{Si)  +  d(Si)<l>(k]Si)<  0. 


(39) 


It  is  possible,  though  tedious  and  not  very  informative,  to  translate  the  conditions  of  The¬ 
orem  3  into  more  explicit  conditions  for  the  linear  stabilizability  of  (1).  This  would  involve 
explicit  limits  on  the  allowed  norm  of  the  control  gain  vector  under  which  the  nominally 
stable  modes  would  be  guaranteed  to  remain  stable  for  the  closed-loop  system.  Obtaining 
such  upper  bounds  from,  say,  Liapunov  matrix  equations  would  be  feasible  but  would  lead  to 
conservative  results.  Given  such  bounds,  one  would  simply  check  if  any  of  the  gain  vectors  of 
Theorem  3  which  result  in  appropriate  values  of  the  stability  coefficients  satisfy  the  bounds. 
We  do  not  present  calculations  on  this  issue  here,  since  they  are  intricate  and  are  not  of 
primary  importance  in  this  work.  We  note,  however,  that  one  approach  to  obtaining  the 
bounds  is  to  change  coordinates  of  the  linearized  system  so  as  to  separate  the  uncontrollable 
critical  mode  from  the  remaining  modes.  Also,  note  that  it  may  be  best  to  obtain  such 
bounds  using  extensive  numerical  search. 

4.3  Construction  of  Stabilizing  Linear  Feedbacks 

Next,  we  combine  the  results  of  the  previous  two  subsections  to  result  in  a  procedure  for  the 
construction  of  a  stabilizing  linear  feedback  u  -  kx  which  achieves  a  prescribed  value  0%  (k)  — 
fd.  For  brevity,  only  case  (a)  of  the  statement  of  Theorem  2  is  addressed.  Construction 
procedures  for  other  cases  follow  similarly. 

Our  strategy  for  constructing  stabilizing  linear  feedback  is  as  follows.  First,  the  family 
of  linear  feedback  gain  vectors  k  achieving  a  prescribed  value  assignment  (d%  =  fd  <  0  (with 
fd{s  =  0  fulfilled)  is  generated  via  a  set  of  linear  algebraic  equations.  Then,  those  k  that 
preserve  the  stability  of  the  open-loop  stable  modes  can  be  selected  from  the  obtained  family. 
If  for  the  current  choice  of  jd  <  0  there  exists  no  such  solution  k ,  decrease  the  modulus  \jd\ 
and  restart  the  procedure. 

Given  a  desired  value  fd  <  0  for  fd%  ,  the  problem  is  equivalent  to  finding  those  k  which 
solve  two  equations:  S(k)  =  Si  and 

=  (40) 

(Note  that  d(Si)  7^  0,  do  to  the  fact  that  A  >  0,  which  holds  because  case  (a)  of  Theorem  2 
is  in  force.) 
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Applying  Lemma  B.l,  we  obtain  the  linear  system 

Ms  =  [A,;,  05,,/?] 


(41) 


where  the  coefficient  matrix  'S' 5  is  given  by 

*5  =  ^s(St,0)  :=  [r  -  SiA-b,A~{B*(r*y)  -  (42) 

Note  that  rank  0)]  =  2  for  any  <St-  and  0  such  that  B*(r*,r*)  7^  <j>s,,pb,  since  Ir  =  1 

whereas  lA~x  =  0  for  any  x  G  Es .  Also,  note  that  if  it  happens  that  P$(Si)  <  0,  then  (39) 
holds  with  (f>{k\Si )  =  0,  giving  a  linear  equation  to  be  solved  along  with  (31). 

Procedure  S  below  summarizes  the  main  steps  in  constructing  a  stabilizing  linear  feedback 
under  hypothesis  (S)  assuming  the  situation  specified  in  case  (a)  of  the  statement  of  Theorem 
2. 

Procedure  S  (Construction  of  stabilizing  linear  feedbacks  in  Case  (S)) 

Step  1.  Solve  Eq.  (31)  for  <Si  and  <S2,  and  compute  rq,  v2,  wi  1  w2 ,  w3,  B*(r*,r*),  P<?(c\) 
and  d(Si)  for  i  =  1  and  2.  Select  any  /3  <  0. 

Step  2.  Compute  4>st,p  for  ?  =  1  and  2. 

Step  3.  Solve  Eq.  (41)  for  k,  first  for  i  =  1  then  for  i  =  2.  If  one  of  these  solutions  k 
satisfies  \\k\\  <  k ,  then  this  k  gives  a  stabilizing  linear  feedback  u  =  kx.  Otherwise, 
select  another  0  <  0  of  smaller  magnitude,  and  go  to  Step  2.  ■ 


5  Linear  Feedback  Stabilization  in  the  Case  of 
One  Pair  of  Imaginary  Eigenvalues 

In  this  section,  results  on  linear  feedback  stabilizability  of  the  origin  of  (1)  in  Case  (H)  are 
given.  The  development  entails  an  analysis  of  the  assignability  of  and  implications  for 
linear  stabilizability.  A  procedure  for  the  synthesis  of  stabilizing  linear  feedback  laws  which 
achieve  a  prescribed  value  0  <  0  of  0^  is  also  given.  The  analysis  proceeds  in  a  fashion 
similar  to  that  for  Case  (S)  in  Section  4.  However,  in  Case  (H)  there  is  no  equality  condition 
analogous  to  the  condition  0*  =  0  in  Case  (S). 


5.1  Stability  Coefficient  Assignability 

A  formula  giving  the  closed-loop  stability  coefficient  0%H  in  terms  of  the  linear  state  feedback 
vector  k  is  obtained  in  Appendix  A.  This  formula  is: 

01,  =  Ph(H)  +  rfi(W)&(*;  W)  +  R e{d2{H)<h(bH)}  (43) 

with  7i  —  7i(k)  the  complex  fraction  given  by  Eq.  (27),  PhOP)  a  real-valued  cubic  polyno¬ 
mial  in  d]  (77)  a  real  coefficient,  c/2(77)  a  complex  coefficient,  <j)\(k]  7i)  a  real  fraction,  and 
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<t>2 (k',Pi)  a  complex  fraction.  These  are  given  by 


Ph{H) 
di(H) 
d2(H) 
<t>  \{k\Pi) 
<f>2(k;7i) 


=  (32H+2Re{H{lwl)  +  H{lw2)  +  H2{lw3)  +  HH(lw4)  +  H2'H(lw5)},  (44) 


=  Re{(lv<n)  +  H(lvn)},  (45) 

=  (^02)  +  Pi(lv  12),  (46) 

=  T{k\A-xB*{r*^\A-xb),  and  (47) 

=  f(k-,(A-2iujcI)-lB*{r*,r*),(A-2iujcI)-lb),  (48) 


respectively;  See  Appendix  A  for  the  values  of  coefficient  vectors  Wi ,  i  =  1, . . . ,  5  and  u0i,  U02, 
voii  vi2  as  well  as  the  vectors  B*(r*,r*)  and  B*(r*,r*).  Note  that  although  Pi  does  not  appear 
explicitly  in  the  formulae  (47)  and  (48)  for  (f>4(k\Pi)  and  cj)2(k]Pi),  it  affects  these  quantities 
through  the  expressions  for  5*(r*,r*)  and  i?*(r*,r*),  respectively.  Compare  with  Eq.  (38) 
for  Case  (S). 

Apparently,  the  fraction  Pi  plays  a  role  similar  to  that  of  S  in  Case  (S),  and  analogies  occur 
between  Ph  and  P$ ;  between  (di(Pi),  d2(Pi))  and  <7(5);  and  between  (<j>i(k]  Pi),  <f>2(k ;  Pi))  and 
<f>(k;S).  However,  because  of  the  absence  of  a  necessary  stabilizability  condition  analogous 
to  /?j  =  0,  Pi  is  a  free  variable.  In  Case  (S),  the  variable  S  was  subject  to  the  requirement 
S  —  Si.  The  difference  between  Eq.  (35)  for  fi2s(k)  and  Eq.  (43)  for  P2fi(k)  results  in  added 
freedom  in  achieving  value  assignments  [P1h  =  fi.  Indeed,  since  Pi  need  not  be  fixed,  /3%h  can 
be  modified  via  three  ‘channels,’  namely  Pi ,  (j>i{k\Pi )  and  (f>2(k\Pi).  In  parallel  to  Theorems 
3  and  4,  we  have 


Theorem  4.  (02h  Assignability)  Let  (1)  satisfy  hypothesis  (H)  and  suppose  lb  =  0.  The 
function  P2h(-)  is  onto  if  i  =  1,2, 3, 5,  as  well  as  /u0i,  lv 02,  lv n  and  lv42  do  not  vanish 
simultaneously.  If  lw{  =  0  for  i  =  1,2, 3, 5  and  lv 0i  =  lvu  =  lv 02  =  ivi2  =  0,  then 

=  (32h  +  |77(^)|2Re{/«;4}  (49) 

which  has  a  global  maximum  (resp.,  global  minimum)  /72j?  if  Rej/u^}  <  0  (resp.,  if  Re{/?<q}  > 

0).  ■ 


Theorem  5  below  is  analogous  to  Theorem  3  of  Section  4.  The  theorem  statement  uses 
the  notation  k  introduced  at  the  end  of  Section  3,  the  notation  1 Z~k(F(-))  introduced  at  the 
end  of  Section  1,  the  notation  and  the  quantities  Pi{k ),  <f>i(k]Pl)  and  <j>2(k;Pi )  given  by  (27), 
(47)  and  (48),  respectively. 

Theorem  5.  (Sufficient  Condition  for  Linear  Stabilizability  in  Case  (H))  Let  (1)  satisfy 
hypothesis  (H)  and  suppose  lb  =  0.  Then  the  origin  of  Eq.  (1)  is  stabilizable  by  linear  state 
feedback  if  there  is  a  k  with  ||L’||  <  k  that  satisfies  kA~1b  ^  —1  and  k(A  —  2iucI)~1b  7^  — 1; 
and  with  Pi  :=  Pi{k ),  we  have  Pi  €  TZ-k(Pi(-))  and  that  <f>i(k;Pi)  €  'TZ-k{<f)\ (•))  and  (j>2{k\Pi)  € 
TZj t(<^2(’))  satisfy  the  inequality 

Ph(PC)  +  d1(Pi)<j)1(k-,  Pi)  +  R e{d,2(Pi)<j>2(k-Pi)}  <  0.  (50) 
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5.2  Construction  of  Stabilizing  Linear  Feedbacks 

In  this  subsection,  which  parallels  Section  4.3,  a  procedure  is  given  for  generating  stabilizing 
linear  feedback  laws  u(x )  =  kx  for  (1). 

An  equation  similar  to  (41)  in  Case  (S)  is  sought  to  characterize  the  set  of  vectors  k 
which  achieve  the  value  assignment  /?£  =  j3.  For  brevity,  as  in  Case  (S),  only  the  simplest 
(nondegenerate)  case  is  considered,  namely  that  in  which  di(J~L)  and  d2(7i)  are  nonvanishing 
functions  of  7i.  That  is,  we  assume  that  at  least  one  of  the  coefficients  appearing  in  each  of 
Eqs.  (45)  and  (46)  does  not  vanish.  Consider  the  assignment  equation 

PH(n(k))  +  d1(H(k))(l>1(k;n(k))  +  R e{d2(n(k))(f>2(k;  H{k))}  =  f3.  (51) 

By  the  reasoning  preceding  Theorem  4,  first  we  hx  an  7i  €  C  such  that  \dl(7{)\  +  \d2(7i)\  / 
0  so  that  @2H  can  be  modified  also  by  either  <f>i(k;H(k))  or  <j)2(k-,'H{k)).  Thus,  Ph{'H ) 
becomes  a  fixed  constant,  and  <f>i(k\H(k ))  and  <j>2(k ,Ti(k))  reduce  to  <f>x(k;H )  and  <f>2(k;H), 
respectively.  Now,  Eq.  (51)  can  be  treated  as  a  linear  equation  in  the  two  unknowns 
(f)\  :=  (f>i{k]  H)  €  K  and  cj)2  :=  <f>2(k;H)  €  C: 

<k{H)<f> r  +  R e{d2(H)<h}  =  0-  Ph{H).  (52) 

It  follows  that  (51)  now  consists  of  three  equations  in  k ,  namely  7i(k)  —  7i,  (f>i(k-,7i)  —  (j>i 
and  (f>2(k]H )  =  (f>2  (with  {<f>\^2)  any  solution  pair  of  Eq.  (52)),  each  involving  expressions 
of  the  form  of  the  fraction  introduced  in  (24).  Finally,  we  use  Lemma  B.2  (including  the 
(-)'-notation  defined  in  Eq.  (B.4)  for  complex  fractions  7Y(A:)  —  7i  and  <f>2 (k;7i)  =  < j>2 )  to 
obtain  the  system  of  real  linear  equations  (see  Eq.  ( B . 5 ) ) : 

k^H  =  [Re{H},  Im{'H};  <j>x;  Re{<^2},  Imj^}]  (53) 

where  the  coefficient  matrix  $ //  is  given  by 

1A2) 

:=  [(r)',  ((A  -  iucI)-b)'-  f*)  -  ^6); 

{{A  —  2iucI)~1  B*(r* ((A  —  2iucI)~1b)'].  (54) 


Procedure  H  below  is  analogous  to  Procedure  S  given  at  the  end  of  Section  4.  It  sum¬ 
marizes  the  main  steps  in  constructing  a  stabilizing  linear  feedback  under  hypothesis  (H). 
The  simplest  situation  addressed  in  the  statement  of  Theorem  4  is  assumed,  namely  that  in 
which  not  all  of  the  quantities  *  =  1, 2,  3,  5  and  /i>oi,  h;o2,  lv n,  lv\2  vanish  simultaneously. 

Procedure  H  (Construction  of  stabilizing  linear  feedbacks  in  Case  (H)) 

Step  1.  Compute  u01 ,  v02,  v, ,,  v12  and  lwl  for  i  =  1, ...  ,5.  Select  any  0  <  0. 

Step  2.  Select  7i  €  C  such  that  dx{7i)  and d2('H)  do  not  both  vanish.  Compute  Ph('H) 
using  Eq.  (44).  Select  a  solution  pair  4>i  6  1R-,  4>2  6  C  to  Eq.  (52). 

Step  3.  Solve  Eq.  (53)  for  k.  If  a  solution  k  satisfies  ||A;||  <  A:,  then  this  k  gives  a  stabilizing 
linear  feedback  u  =  kx.  Otherwise,  select  another  /A  <  0  of  smaller  magnitude,  and  go 
to  Step  2.  ■ 


6  Conclusions 


Linear  feedback  stabilizability  and  stabilization  of  nonlinear  systems  with  uncontrollable 
critical  modes  have  been  studied.  The  paper  focused  on  the  two  basic  critical  cases,  those 
in  which  the  system  linearization  at  the  origin  possesses  either  a  simple  zero  eigenvalue  or 
a  pair  of  simple  pure  imaginary  eigenvalues.  Construction  procedures  facilitating  synthesis 
of  stabilizing  linear  feedback  laws  were  also  given.  In  general,  the  conditions  obtained  com¬ 
bine  sufficient  conditions  for  assignability  of  certain  associated  stability  coefficients  with  a 
requirement  that  open-loop  stable  eigenvalues  not  be  destabilized.  The  critical  system  is 
generically  stabilizable  by  a  linear  static  state  feedback  u  =  kx ,  under  the  condition  that 
the  open-loop  stable  eigenvalues  have  a  sufficiently  large  margin  of  stability.  In  the  case 
of  systems  (1)  for  which  the  control  u  enters  the  vector  field  f(x,u)  linearly,  the  results 
presented  here  readily  simplify,  since  destabilization  of  stable  modes  is  no  longer  a  possi¬ 
bility.  The  specialized  results  for  these  so-called  linear-analytic  systems  were  not  given  in 
this  paper,  in  the  interest  of  brevity.  An  interesting  topic  for  further  research  concerns  the 
trade-off  between  between  linear  and  nonlinear  feedback  in  the  local  stabilization  of  systems 
for  which  the  critical  eigenvalues  are  nearly  uncontrollable. 

Appendix  A.  Closed-Loop  System  Stability  Coefficients 

Here,  we  compute  the  stability  coefficients  for  the  closed-loop  system  (13)  under  hypotheses 
(S)  and  (Id).  These  are  /3js  and  in  Case  (S),  and  in  Case  (II).  To  employ  the 
general  expressions  for  stability  coefficients  given  in  Facts  1  and  2,  it  is  necessary  that  the 
bilinear  and  trilinear  forms  in  the  closed-loop  dynamics  (13)  retain  the  symmetry  properties 
of  the  analogous  terms  in  the  open-loop  dynamics  (4).  To  achieve  this,  we  employ  the 
symmetrization  operation,  as  in  (Abed  and  Fu,  1986),  (Abed  and  Fu,  1987). 

For  the  quadratic  and  cubic  forms  B*(x ,  x )  and  C*(x ,  x,  x),  respectively,  of  Eqs.  (15)  and 
(16),  the  associated  bilinear  and  trilinear  forms  are  defined  as  follows: 

B*(x,  y)  :=  B(x,y)  +  \[{kx)Anj  +  (ky)A1x]  +  b2(kx)(ky ),  (A.l) 

C*(x,y,z)  :=  C(x,y,z)  +  J[(Ar®)5i(y, z)  +  (ky)B1(z,x)  +  ( kz)B1(x,y )] 

-\-\[{kx)(ky)A2z  +  (ky)(kz)A2x  +  ( kz)(kx)A2y }  +  b3(kx)(ky)(kz).  (A.  2) 

We  can  now  proceed  to  the  calculations  associated  with  obtaining  the  expressions  for  the 
stability  coefficients  ,  (3% .  occurring  in  the  discussion  of  Case  (S).  To  begin,  we  compute 
the  coefficient  vectors  iq  and  v2  first  appearing  in  the  text  in  Eq.  (30).  Letting  x  :=  r*,  with 
r*  as  given  in  (25)  (Case  (S)),  in  Eq.  (15)  yields  Eq.  (30),  with  the  coefficient  vectors  iq 
and  v2  given  by 

:=  -2B(r,A~b)  +  Apr,  (A.3) 

v2  :=  B{A~b,A~b)  -  AiA-b  +  h,  (A.4) 

respectively.  Now  using  this  same  substitution  x  :=  r*  in  Eq.  (16)  yields  C*(r*,  r*,  r*), 
a  quantity  required  in  applying  Eq.  (17b).  The  vector  A*~ B*(r* ,  r*)  which  occurs  in  Eq. 
(17b)  is  easily  obtained  using  Eqs.  (22)  and  (30).  Substituting  this  value  and  that  of  r* 
into  (A.l)  gives  the  first  term  in  (17b).  Thus,  (3*ls  is  given  by  (35),  along  with  the  notation 
(36)-(38),  where 

w0  :=  -2 B(r,A~B(r,r))  +  C(r,r,r), 

uq  :=  Bi(r,r )  —  A\A~ B(r,r)  —  2 B(r,A~b)  —  3 C(r,r,  A~b), 
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(A- 5) 
(A- 6) 


w2  :=  A2r  —  2Bx(r,A  B)  +  3C(r,r,A  b),  and 
w3  :=  b3-C(A~b,A-b,A-b). 


(A.7) 

(A.S) 


Next  we  obtain  the  closed-loop  stability  coefficient  /32h  which  is  defined  in  Case  (H). 
First,  using  r*  from  Eq.  (25)  (Case  (H)),  substitute  x  r*  and  y  f*  into  Eq.  (A.l),  and 
substitute  x  :=  r*  into  (15),  to  obtain 


B*(r*,  f*)  =  B(r,  r )  A  Tiv\  A  7fwx  A  'H'Hv2, 

(A. 9) 

B*(r*,  r*)  =  B(r ,  r)  A  Hv 3  A  H2vA. 

(A.10) 

To  state  the  values  of  the  vectors  vx,v2,v3,V4  appearing  in 

(A. 9)  and  (A. 10),  denote 

£  :=  (A  -  iucI)~b. 

(A.ll) 

Then  the  u,-,  i  —  1, . . . ,  4  are  given  by 

ux  :=  -5(r,£)  A  \Axf , 

(A.12) 

v2  :=  B(£,  £)  —  ^Aj(£  A  £)  +  b2, 

(A- 13) 

v3  :=  —2 B(r,  £)  A  Axr,  and 

(A- 14) 

u4  :=  -B(£,£)  —  Ax£  A  b2. 

(A. 15) 

Next,  substitute  x  :=  r*,  y  :=  r*  and  z  :=  r*  into  Eq.  (A. 2)  to  obtain  C*(r*,  r* ,  f*).  Then, 
premultiply  the  quantities  B*(r*,  r*)  and  B*(r*,  r*)  by  A*-1  and  (A*  -  2?'uy)~1,  respectively. 
(These  matrix  inverses  are  given  in  Eqs.  (28)  and  (29),  respectively.)  Substitute  the  obtained 
products  along  with  r*  and  r*  in  place  of  the  proper  vector  arguments  in  (A.l)  to  obtain  the 
first  two  terms  in  the  expression  (18)  for  /32}j.  The  final  expression  (43)  follows  by  adding 
these  two  vectors  to  C*(r*,  r*,  r*),  obtained  earlier.  The  notation  used  in  Eqs.  (44)-(48) 


(needed  in  the  expression  (43))  is  as  follows: 

t>oi  :=  -2 B(r,  A~lb)  A  Ayr,  (A. 16) 

vn  :=  2 B(z,  A~lb)  —  Ai(z  +  A~lb )  +  2 />2,  (A.17) 

v02  :=  -B(f,  A-1!,)  +  iAjr,  (A.18) 

V12  :=  B(z,  A-1!))  -\Ax{z  +  A-1 5)  +  b2,  (A.19) 

wx  :=  —  B(r,  A_1rq)  —  2 B(z,  B{r,r))  —  \ B(r,(A  —  2iu>cI)~'lv2 ) 

-f  C(r,r,z)  A  Am  A  §£x(r,r),  (A.20) 

w2  :=  — B(r,  A_1ui)  —  B(z,b) 

-Jtf(r, r,z)  A  JA,6  A  Jj?i(r,r),  (A.21) 

u>3  :=  5(z,  A_1ux)  —  l/?(f,  (A  —  2fu;c/)_1i;4) 

+\C(f,z,z )  -  lAiA_1ui  -  \Bx(r,z)  +  \A2f ,  (A. 22) 

W4  :=  — S(r,  A_1u2)  +  |i?(^,  A~lvx)  A  |C(r, s) 

—  ^AxA-1';)!  —  \AX(A  -  2iu>cI)~1v2  —  |F?x(r,  z)  +  |A2r,  and  (A. 23) 

w5  :=  B(z ,  A_1u2)  +  |f7(z,  (A  —  2jo>c)-1u4) 

+fC(.z,  z,  z)  —  |AXA  ju2  —  ^AX(A  —  2iucI)  1u4 

+~BX (z,  z  A  2 z)  - \A2(2z  A  z)  A  f  k-  (A.24) 


Appendix  B.  Solution  Sets  of  Certain  Rational  Equations 

In  this  Appendix,  we  consider  the  solution  of  systems  of  equations  in  which  certain  rational 
functions  of  the  unknown  vectors  appear.  Each  rational  function  is  a  ratio  of  linear  functions 
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of  the  unknown  vectors.  The  systems  of  equations  are  related  to  systems  of  linear  algebraic 
equations,  and  the  structure  of  the  solution  set  is  studied. 

Consider  the  following  special  class  of  functions  F  given  by  the  rational  functions  of  a 
real  vector  k ,  which  are  parametrized  by  real  vectors  u  and  v: 

^(M,>0:=TTk'  <BI> 

The  range  of  F(k]u,v),  and  the  set  of  solutions  k  to  a  system  J-{k]u,v)  —  9  (9  £  IR),  are 
important  in  the  stabilizability  analysis  and  stabilization  procedures  in  the  paper. 

First,  we  note  that  for  any  9  the  equation  F(k-,u,v)  =  9  is  equivalent  to  the  linear 
equation  k(u—9v )  =  9.  Hence,  for  any  6  such  that  u—Ov  ^  0,  the  solution  set  of  T{k\  u,  v)  =  6 
coincides  with  that  of  the  linear  equation  k(u  —  9v)  =  9.  (Note:  thus,  if  u,v  are  linearly 
independent,  the  equation  Jr(k',u,v)  —  9  is  solvable  for  any  9 ;  if  u  =  av  for  some  a  9, 
then  u  —  9v  =  (a  —  6)v  7^  0  and  the  rational  equation  is  solvable  for  any  9  ^  a.)  Applying 
this  idea  to  a  system  of  equations  each  in  the  form  of  (B.l):  J-t(k\  u,  v )  =  0,-,  i  —  1, . . . ,  n, 
we  have 

Lemma  B.l.  The  solution  set  of  the  simultaneous  equations  J~i(k]  w,-,  u,)  =  9i,  i  =  1, . . . ,  n 
defined  in  the  form  of  (B.l)  via  a  set  of  vectors  {ui,  iq, . . . ,  um,  vm},  coincides  with  the 
solution  set  of  the  linear  system 

&[t/i  9\V\, . . . ,  um  -  9ni cm ]  —  [^i)  •  •  •  (B.2) 


Note  that  if  kv{  =  —1  for  some  solution  k  of  (B.2)  for  some  i ,  then  kul  =  0  necessarily. 
Also  note  that  (B.2)  is  guaranteed  solvable  for  any  $i, . . . ,  if  the  set  {(q  —  9\V\, . . . ,  um  — 
is  linearly  independent. 


The  conclusion  above,  in  particular  whether  or  not  a  single  rational  function  of  type 
JF(-;  u,v)  is  onto,  is  in  doubt  if  an  additional  constraint  is  imposed  on  k ,  such  as  |&|  <  k. 


Direct  application  of  Lemma  B.l  in  Case  (H)  is  not  possible,  due  to  the  presence  of 
complex  quantities.  Therefore  some  modification  in  the  statement  of  Lemma  B.l  is  in  order. 
Consider  Eq.  (B.l),  with  the  function  F(k;u,v)  now  defined  for  vectors  u,v  6  C".  For 
brevity,  denote  by  [-]r  and  [-]1  the  real  and  imaginary  parts  of  a  vector  [•].  Rewrite  A  ( k. ;  «,  w) 
as: 


F(k]  'u ,  v) 


( kuT )  +  i(ku’) 

(1  +  kvr)  +  i(kv1) 


(B.3) 


Let 

I  («)'  :=  1,1  ~  Re{fl}er  +  Imj^}?;1,  ,R  ., 

I  (v)'  :=  v'  -  1  ’ 

It  can  be  shown  that  for  any  9  £  €  for  which  the  set  {(«)',  (u)'}  is  linearly  independent, 
or  lm{0}(u)'  =  Re{0}(u)',  the  equation  F(k;  u,  v)  —  9  has  solutions  coinciding  with  the 
solutions  of  the  system  of  linear  equations  k[(u)\  (u)']  =  [Re{#},  lm{0}]. 
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The  same  idea  applies  to  a  set  of  rational  equations  lFi(k;iii,Vi )  =  0,  , i  =  1  < 

n/2,  all  in  the  form  of  (B.l)  defined  via  vectors  That  is,  to  consider  each  equation 

ui,Vi)  =  0{  as  in  (B.3)  and  to  define  {uf)'  and  (v,-)7  as  in  (B.4).  We  thus  have 

Lemma  B.2.  The  solution  set  of  the  rational  system  T{{k\  Ui ,  i>i)  =  6i  £  C,  i  =  1, . . . ,  m  < 
n/2  coincides  with  the  solution  set  of  the  linear  system 

^[(wi)',  (wi)';  •  •  ■ ;  ( Um)' ,  («m)']  =  [Re{6>j},  Im-fffi}; . . . ;  Re{6>m},  lm{0.m}].  (B.5) 
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